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Abstract. We study random Hamiltonians on finite-size cubes and waveguide segments 
of increasing diameter. The number of random parameters determining the operator is 
proportional to the volume of the cube. In the asymptotic regime where the cube size, and 
consequently the number of parameters as well, tends to infinity, we derive deterministic and 
probabilistic variational bounds on the lowest eigenvalue, i. e. the spectral minimum, as well 
as exponential off-diagonal decay of the Green function at energies above, but close to the 
overall spectral bottom. 


1. Introduction 

Quantum disordered systems often exhibit localization, i.e. the absence of propagation of 
wavepackets. For random ergodic Schrddinger operators in L2(R"') this has been established 
in various regions in the energy x disorder diagram. For such models, localization comes about 
thanks to the local effect of random variables (encoding the disorder in the Hamiltonian) and a 
global conspiracy of randomness over large scales. A natural approach to study, and actually, 
prove localization, is to analyze first the spectral effects of a single random variable on a 
specific type of random operator, then the cumulative effect of many variables on finite but 
large cubes in configuration space, and finally conclude that a quantitative form of localization 
persists if one takes the macroscopic limit. 

We take a reverse, conceptual and abstract approach. We want to formulate criteria on 
the properties of local perturbations (single site potentials for usual random Schrddinger 
operators) which ensure that localization will ensue in an appropriate disorder/energy regime. 
To illustrate what we mean, let us consider the very first result on localization in L2(R”') 
obtained by Holden and Martinelli in m- They consider the random Schrddinger operator 

= — A + '/Z ujk u{- — k) in L2(R"'), where ojk, k € 'Iff, are uniformly distributed on [ 0 , 1 ] 
fceZ'i 
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and u{x) = X[o,i]"(3^) is the characteristic function of the unit cube. The global strategy 
employed in m to prove localization is the multiscale analysis of Frohlich and Spencer. On 
the local level, the properties of the single site perturbation oj^ eA ujk u{- — k) are essential. It 
is linear in nonnegative on L 2 (R"') and strictly positive on L2([0,1]”). The question we 
raise is: If the function u = X[o,i]'i is replaced by a more general function, or even an operator 
distinct from an multiplication operator, which properties should it have in order to ensure 
pure point spectrum of While this has been studied for a number of specific, physically 
relevant, models, our approach is conceptual. We want to understand a set of sufficient 
conditions on the local building blocks of the Hamiltonian (single site perturbations) which 
ensures localization, or at least important partial results used on the road to localization. 
Although in this paper we do not provide a complete answer to the above question, we make 
a first important step. Namely, we provide an initial length scale estimate, that is one of the 
main steps in proving spectral localization via multiscale analysis, for a very wide class of 
random Hamiltonians with weak disorder. 

The second key ingredient to make the multiscale analysis work is a Wegner estimate. The 
role of the two ingredients is the following: the multiscale analysis is a induction procedure 
over a sequence of increasing length scales. While the initial length scale estimate provides 
induction anchor, the Wegner estimate guarantees that the induction step works. Physically, 
the Wegner estimate ensures that resonances between spectra of disjoint subsystems occur 
only with small probability. In a sequel paper we plan to give a set of conditions on general, 
abstract random Hamiltonians which imply the Wegner estimate. This set is distinct, but 
similar to the conditions we impose in the present paper. Thus for random Hamiltonians 
which satisfy both requirements localization via multiscale analysis follows. 

The indication how to implement the proof of Wegner estimate is provided by Lemma 
2.3 below. It describes the lifting of the spectral bottom for periodic configurations of the 
random coupling constants. This ensures that there is a (small) energy interval near the 
minimum of the spectrum of the original, unperturbed operator which is uncovered by the 
random perturbations: There exists a operator in the ensemble whose resolvent set contains 
the mentioned energy interval. In this situation the vector-field method introduced in 

[T6] by Klopp and developed in m and m can be applied. 

While our theorems cover a substantially more general setting, let us describe here in the 
introduction a special case of the model we consider: Let Ci,C 2 ■ 1]"") —)• L2([0,1]"") be 

bounded symmetric linear operators, e > 0, C{t) := tCi + t^C 2 , t G [—e,e], be an operator 
family, {S{k)u){y) = u{y + k),y G H”, k £ he the shift operator, k G Z”, be a sequence 
of numbers with values in [—1,1], and 

:= —A -I- := ^^S{k)C{£UJk)S{—k), 

fcer 

We have no specific requirements on the type of operators Ci,C 2 - they could be multiplication, 
differential, or integral operators or a combination of these. In Section [3| we give a number of 
examples covered by our general setting, including scalar potentials, magnetic fields, random 
metrics, Laplacians on infinite strips and layers with random boundary, as well as integral 
operators. 

Our two assumptions on the single site operators C,\^C ,2 are the following: Let —A be the 
negative Neumann Laplacian on [0,1]"', 1: [0, Ij” —)• 1 the constant function and u the unique 
solution to — Au = £il with f udy = 0. We assume that 

[ 0 , 1 ]" 
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(Al’) f Ci\dy = 0 and (A2’) f C 2 ^dy — f uC]\dy > 0. 

[ 0 , 1 ]" [ 0 , 1 ]" [ 0 , 1 ]" 

Our first result is that the lowest eigenvalue X%{co) of the restriction of to 

n 

nAr := {y € m"" ■ y = S (0, A^)} with Neumann boundary conditions obeys 

i=\ 

X%{u:) ^ C 2 e^ Y. := Z" n [0, N)^ 

fceTjv 

provided N > Ni, and 0 < e < ^. Here ci, C 2 , A^i S (0,oo) are independent of e and N. 

If cjfc, k £ form an i.i.d. sequence of random variables, we deduce that for small, but not 
too small values of e > 0 

where constant C 4 > 0 depends only on the distribution of loq. Finally, we prove a Combes- 
Thomas bound for the general class of operators we consider, and derive a initial scale esti¬ 
mate, as it is used for the induction anchor of the multiscale analysis. A Combes-Thomas esti¬ 
mate is a bound on the off-diagonal exponential decay of the integral kernel of — E)~^ 

provided an a-priori lower bound on dist(i]l, a{J-L%{uj))) is known. Here cr(-) denotes the spec¬ 
trum of an operator. The novelty of our result is that the considered operators need only be 

block diagonal with respect to the decomposition 0 T 2 ([ 0 , 1 ]"" + k'), but not necessarily a 

fceZ" 

differential operator. 

In Section [3.71 we show that the model described in this introduction is covered by the more 
general, abstract model defined in Section [2j 

History and earlier results. 

Results on random waveguides. The results presented here are a generalization and improve¬ 
ment of those in [3] : In [3] we considered randomly wiggled quantum waveguides in the ambient 
space IR^. For this specific model a variational estimate analogous to Theorem 12.II was estab¬ 
lished in [31 Corollary 4.2]. In [3] we studied an apparently very similar disordered model, 
namely a randomly curved waveguide. For this model the hypotheses (Al) and (A2) are not 
satisfied and our analysis showed that the lowest eigenvalue exhibits a behaviour dis¬ 

tinct (in some sense opposite) to the one encoded in inequality ()4.2|] . Common to both types 
of random waveguides studied in |3] and |3] is a non-monotone dependence on the random 
variables Uk- This is a challenge to the mathematical analysis, as will be elaborated further 
below. To the best of our knowledge the first disordered model of a quantum waveguide was 
studied in m- There the width of the waveguide is determined by a sequence of random 
parameters Uk- This gives rise to a monotone influence of the parameters and facilitates the 
study of the spectrum. In m in addition to an initial length scale estimate for the Green’s 
function a Wegner estimate was provided, yielding spectral localization near the bottom of 
the spectrum. 

Weak disorder quantum Hamiltonians. Our model depends on a global parameter e > 0. It 
tunes the overall strength of the disorder present in the operator. The interest is now, to 
identify an energy interval, depending on the disorder strength e, where an initial length 
scale estimate, a Wegner estimate and spectral localization hold. Corollary 12.71 provides such 
a statement concerning the initial length scale estimate. Most detailed results identifying 
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energy regimes with spectral localization in the weak disorder regime have been obtained for 
the Anderson model on or its continuum analog, the alloy type model on L^(]R"'). 

Corresponding to the general setting of this paper, we will restrict our discussion to continuum 
models, i.e. quantum Hamiltonians dehned on H” or an open, unbounded subset thereof. In 
|18] it is proven that under the assumption 

(-A1) u G L^(IR”), J uix) dx^O 

there is an energy interval G IR, with size of the order e such that the alloy type model 

— A + V^er + e 'w(' ~ k) 

fceZ" 

exhibits spectral and dynamical localization in almost surely. Here Vper is a bounded Z”- 
periodic potential. Due to assumption (-lAl), [18] does not cover the situation considered 
here (if we assume that the random perturbation is a multiplication operator). In this sense 
our result, when specialised to the case that the random part of the Hamiltonian is a potential, 
complements the result of |18j . 

Spectral analysis of non-monotone random Hamiltonians. The proofs of initial length scale 
estimates and Wegner estimates simplify greatly if the random variables cok influence the qua¬ 
dratic form associated to the Hamiltonian in a monotone way. If this monotonicity property 
is violated one has to identify and use specihc properties of the model at hand in order to 
replace monotonicity. This has been carried out for alloy type models of changing sign e.g. in 
[ITlESKISlESKislEaEniEel, for random displacement models e.g. in [EKDEailHI, for 
random magnetic fields e.g. in ISSlEolliaEIlEIlEI, PE], and Laplace-Beltrami operators 
with random metrics e.g. in |25( [23 ] I24j. 

Innovations. We list the results, methods and conceptual innovations obtained in the paper. 

• We establish a variational lower bound for the spectral minimum of random Hamil¬ 

tonians an arbitrary large, finite boxes H^v. As the box size N grows, the number of 
random variables influencing the random Hamiltonian grows as well, namely 

proportional to the volume of Hv- Thus the variational problem involves a large (and 
increasing) number or parameters. 

• The basic assumption on the influence of the individual parameters on the random 
Hamiltonian T-L^{uj) is the validity of a certain Taylor formula, cf. (12.11) . as well as 
Assumptions (Al’) & (A2’), or their generalizations (Al) &: (A2). In contrast to the 
standard approach, we do not require the dependence to be linear (not even rational). 

• The variational bounds are proven using a perturbative framework based on a non¬ 
self-adjoint modification of Birman-Schwinger principle proposed and developed in 
m, see also 121, El, and |5|. This abstract approach allows a uniform treatment 
of many types of random Hamiltionians studied before (random scalar potentials, 
random magnetic fields, randomly perturbed quantum waveguides), as well as new 
types (e.g. integral operators, randomly perturbed quantum layers). 

• We establish a general Combes-Thomas estimate. It does not require the Hamiltonian 
to be a differential operator, rather it could contain an integral operator part as well. 
To the best of knowledge of the authors such estimates have been so far obtained only 
for local operators. 
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• Now two probabilistic results follow: First we establish an upper bound on the prob¬ 

ability of hnding an eigenvalue of "^^( 0 ;) very close to Aq. The probability is expo¬ 
nentially small in the size N, while the notion of ’very close’ depends on N as well. 
With the same probability we establish that the Greens function for energies above, 
but close to Aq of the resolvent of decays exponentially in space. 

• The size of the energy interval above Aq can be expressed as a function of the weak 
coupling parameter e, instead as of N. This implies an estimate on the ’Lifschitz tail 
regime’ quantified in terms of a small disorder parameter e. The size of the interval 
is not quite, but pretty close to quadratic in s. A quadratic behaviour is the best one 
could expect. 

2. Formulation of problem and main results 

Let x' = {xi,... ,Xn), X = (x',Xn+i) be Cartesian coordinates in IR"' and IR^+^j respectively, 
where n ^ 1. By 11 we denote the multidimensional layer 11 := {x : 0 < x^+i < d} of width 
d > 0. In space IR” we introduce a periodic lattice F with a basis ei, ..., e^; the unit 

n 

cell of this lattice is denoted by i.e., □' := {x' : x' = G (0)1)}- We denote 

i=l 

□ X (0,d). 

For some to > 0 we denote by £(t), t G [—toj^o]) a family of linear operators from 
into ^ 2 ( 0 ) given by 

(2.1) C{t) ■.= tCi+fC2 + t^C^{t), 

where Ci : —)• L^ip) are bounded symmetric linear operators and T 3 (t) is bounded 

uniformly in t G [—to,to]- 

Given u G iL^(n), it is clear that u G and function Ciu is thus well-dehned as an 

element of L 2 (TI). Now we can extend the function CiU by zero in 11 \ □ and this extension is 
an element of L 2 (n). In the sense of the above continuation, in what follows, we regard the 
operators Ci as acting from into L 2 (II). We stress that, in general, the operators Ci 

are unbounded as operators in L 2 (II). 

The main object of our study is the operator 

(2.2) n^{uj) := -A + Bo + T"(w), £"(a;) := ^ S{k)C{£UJk)S{-k), 

fcer 

in n. Here e is a small positive parameter, 0 ;^, k G Z” a sequence of numbers with values in 
[—1,1], Vo(x) = Vi){xn+i) is a measurable bounded potential depending only on the transversal 
variable x^+i, S{k) stands for the shift operator: {S{k)u){x) = u{x' + k, x^+i). The boundary 
condition on clH is either of Dirichlet or Neumann type. We denote this condition by 

(2.3) Hu = 0 

on (9n, and Bu = u or Bu = We consider also the situations when on the upper and 

lower boundaries of SH we have different boundary conditions. Say, on the upper boundary 
we have Dirichlet condition, while on the lower boundary Neumann condition is imposed. 

We consider the operator kC as an unbounded one in T 2 (n) on the domain := {u G 

iL^(n) : (12.31) is satisfied on SH}. The action of the second term in the right hand side of 
(12.2p can be also understood as follows: Given u G we consider the restriction of u on 

the cell Dfc := {x : x — (A:, 0) G □} for each A: G F. Then, identifying cells Dfc and □, we apply 
the operator C{£UJk) to and the result is how C^u is defined on Dfc. 
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For sufficiently small e operator is relatively bounded w.r.t. the Laplacian on T){T-L^) 
with relative bound smaller than one and the latter operator is self-adjoint. Hence, by the 
Kato-Rellich theorem operator T-L^ is self-adjoint, as well. 

Our results concern operators on large, finite pieces 

n 

(2.4) Ho^at := jx : x' = a+ '^ aiCi, a* G (0, N), 0 < Xn+i < d|, 

i=l 

of the layer H, where a G F and G IN are arbitrary. We let 

n 

^a,N := |x' G F : x' = a + ^ ajCj, Oj = 0,1,..., - l| 

i=l 

and observe that Hq^at = U nfe. 

We introduce the operator 

+ ^0 +'^a,Ar(w), C% p^(u):= ^ S{k)C{£tOk)S{—k) 

k&a,N 

in A 2 (nQ,,Ar) subject to boundary condition (12.3p on ^a,N '■= SHa^jv O SH and to Neumann 
condition on ^Hq^at \ The domain of Af(‘^) i® 

(2.5) Ti('H^ Af) ■= G H^(^a,N) '■ u satisfies (12.3h on 7 o, 7 v 

and Neumann condition on SH^^Af \ 7 o,v}- 

The reason why we impose Neumann boundary conditions is the following: We want to give 
lower bounds on the first eigenvalue of finite volume Hamiltionians Af(^)- Since Neumann 
conditions produce the lowest ground state energy, this covers the ‘worst case scenario’. 

By Aq we denote the lowest eigenvalue of the operator 

^- \-Vo on ( 0 ,d) 

“®n+l 

subject to boundary condition (12.311 . The associated eigenfunction normalized in L2{0,d) is 
denoted by ipo : (0, d) —)• IR. Let Tin be the Schrodinger operator —A -|- Vq in T1 subject to 
boundary condition (12.3p on 50 H SH and to Neumann condition on SD \ SH. Note that 
^□V’o = ^oV’O) where here i/’o • Tl —)• R is given by the longitudinally constant extension 
ipoix) = 'ijjo{xn+i)- The second-lowest eigenvalue of Tin is denoted by Ai. 

We make the following assumptions for operators Tj. 

(Al) The identity 

(£iV'o,'</'o)L 2 (n) = 0 

holds true. 

(A2) Let U be the solution to the equation 

(2.6) {Tin - Ao)U = Ci^Po, 
and orthogonal to ipo in L2{0). We assume that 

(2-7) Co := (£ 2 V’ 0 ) V’o)i, 2 (n) “ (t7,TiiAo)L2(n) > 0- 

The two conditions on U in Assumption (IA2p are uniquely solvable since by Assumption (lAip 
CiijjQ is orthogonal to ipo in L 2 (n). 

By K,n we denote the smallest eigenvalue of Ti^ Af- result reads as follows: 
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Theorem 2.1. There exist positive constants ci, C 2 , Ni such that for 

(2.8) N ^ Ni, and 0 < e < 
the estimate 

^SPa.Ar 

holds true. In partieular, the minimum of Xf ^ w.r.t. tOk is Aq and it is aehieved as ujk = 0, 
k G 

Remark 2.2. A very simple interpretation of Assumptions (jAl|) . ()A2p is as follows. They are 
equivalent to the condition 

(2.9) There exists C = const >0 : A* — Aq ^ for t G R, \t\ small . 

for the lowest eigenvalue A* of the operator "Hn + ^(i)- The reason is that the three-term 
asymptotics for A* reads as 

A* = Aq -I- t{Ci'ifo, '0o)i,2(n) + V'o)L2(n) “ {U, TiV’o)L2(n)) + 0{t^). 

Hence Assumptions ()Alh . (IA2I) are equivalent to (12.9p . 

Inequality (12.9p yields that the minimum of A* w.r.t. t is Aq and it is achieved at f = 0. 
In Theorem 14.11 (see also Lemma 14.31) we prove the same for jy, i.e., A^ jy is minimal as 
the perturbation is absent. And this happens mostly thanks to Assumption (jAip . There are 
similar but distinct models, where minimizing the ground state eigenvalue corresponds not 
to the minimal (i.e. absent) perturbation, but to the maximal one, cf. [4j. 

We observe that in order to satisfy Assumption (jA2ll . the scalar product (T2V’0)'0o)L2(n) 
must be positive. The reason is that (ff, Ti'0o)l2(d) ^ 0- Indeed, integrating by parts and 
applying the minimax principle, it is easy to see that 

(t/,AV’o)L2(n) = l|VC/|lL(n) + (^oC/,C/)L2(n) - Ao||C/||i^(n) ^ (Ai - Ao)||C/||i^(n), 
where Ai is the second lowest eigenvalue of TLu, since U is orthogonal to fjQ. This inequality 
provides also an upper bound for ([/,£i^o)i, 2 (n)' First it implies 

II^A||L2(n) < l|Ti'0o||L2(a) 

and thus, 

|(17,£iV'o)L2(n)| ^ _ ^^ll^i^o|li2(n)- 

Then a sufficient condition ensuring (Eip is 

(£2'0o,'0o)L2(n) > ll^i'*/’o|li2(n)- 

Property (12.9p implies an estimate on the spectral minimum of the operator TL^{oj) on the 
infinite domain for periodic configurations w G H. 

Lemma 2.3. Consider the particular configuration w G H with = 1 for all k £ T. Then 
there exists p G (0, oo) independent ofe,a,N, such that 

VQ!Gr,A^GlN,£^0: A^ jy(ti;) ^ Aq -|- cq£^ — pe^ 

and 

Ve^O: inf cj(^^(a))) ^ Aq-I- coe^ —/oe^. 
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Our last deterministic result provides a Combes-Thomas estimate for the general class of 
operators we consider. 


Theorem 2 . 4 . Let a,/3i,/32 G T, mi,m2 G IN 6e such that Bi := C B2 := 

n^2,m2 nc^AT. There exists N2 G ^ such that for N N2 the hound 

(2.10) llxs.(WiN{i^) - A)-AB.IIi,,n„„)^L,(n.,„) < 

holds, where Ci, C2 are positive constants independent of e, a, N, 6 , j 3 i, ( 32 , 1^1, 'm2, ^ CLnd 
5 := dist{X, a{ni j,j{ui))) > 0 . 

Now we formulate our probabilistic results, and introduce for this purpose the assumptions 
on the randomness. Let oj := {cufcjfcgr be a sequence of independent identically distributed 
random variables with the distribution measure p., with support in [— 1 , 1 ]. We assume that 
6_ ^ 0 ^ 6+ and 6_ < where 6_ = minsupp/x and 6+ = maxsupp/x. This gives rise to 
the product probability measure P = fbe configuration space 11 := Xfcgr[~l)l]; 

the elements of this space are sequences u := {cufcjfcgr- By E(-) we denote the expectation 
value of a random variable w.r.t. P. 

Now we are in position to formulate our main probabilistic results. 


Theorem 2 . 5 . Let 7 G IN, 7 ^ 17 . Then for N ^ Ni, where Ni comes from Theorem \ 2 . 1 l 
the interval 


In ■= 


C 3 _ 

E(|a;fc|)iV4 jVG 


C3 := 


'/^2 


is non-empty. For N ^ Ni and e G In, the estimate 


P G H - Ao ^ ^ 

holds true. Here the constant C4 > 0 depends on /x only. 
Our next statement is the initial length scale estimate. 


Theorem 2 . 6 . Let a G T, 7 G IN, 7 ^ 17 , A" ^ W; (in-d e G In- Fix /3i,/32 G Ta^jV; 
mi,m2 > 0 such that Bi := C LIc^jV; B2 := n^2,m2 Then there exists a 

constant C5 independent of e, a, N, / 3 i, /32, mi, m2 such that for N ^ msLx{Nif, Kj, N2} 


V A ^ Ag + 


IIXBi('^a,A - ^) ^XB2 


2y/N 


^ 2\/iVe 


C 5 dist(B]^,B 2 ) 
Vn 


>l-N 


1— 


^-CiN T 


where || • || denotes the norm of an operator in L2{Ila,N) and XB stands for the characteristic 
function of set B. 


Corollary 2 . 7 . Let a G T, 7 G IN, 7 ^ 17 and ci be as in i 2 . 8 \) . Choose e > 0 and A G IN, 
A ^ max{N^, Kf, N2}, such that A = {e/ci)~'^^^. Let /3i,/32, mi, m2, Bi, B2, C5 be as in 
Theorem lil. 61 Then 


P VA ^ Ao + 


1 / e A 


7/8 


) : IIXBi('^a,V - A) 


2 \cij 


-7/8 


-C5 dist(Bi,B2)(^) 


\ 7 / 8 '' 


-n(7-l)/2 


— C4 — 

e 


^ 1 - 


e 

Cl 
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Note that, since N ~ e, and since in applications we have dist(i?i, i? 2 ) ~ iV we have 

/ P \ 7/8 ^ 

dist(5i,52) i—j ~ ~ A^2 > 1 

Thus we are indeed witnessing an off-diagonal exponential decay of the Green’s function, with 
high probability. With the smallest value 7 = 17 which is allowed, we obtain an energy interval 
with width of order which is not much smaller than Intuitively, one would expect 
that in the weak disorder regime (if the first order perturbation annihilates) the Lifschitz tail 
regime interval is of order So, in this respect our result is suboptimal. This is the price 
we pay for treating very general perturbations, instead of, say, just multiplication operators. 

Corollary 12.71 quantifies a Lifschitz-tail regime (an energy interval) in the weak disorder 
regime (a small constant multiplying the random variables). Lifschitz tails denote the ex¬ 
ponential thinness of the infinite volume integrated density of states near the bottom of the 
spectrum. Our results have nothing to say about the integrated density of states, since in the 
limit N ^ 00 the coupling e shrinks to zero. However, when it comes to proving localization, 
one always uses some kind of finite volume criterion, like the probabilistic initial length scale 
decay estimate for the Green’s function (even is Lifschitz asymptotics of the integrated den¬ 
sity of states have been established). Thus, for this purpose our estimate is equally good as 
establishing a Lifschitz tail regime on the energy axis. 

Remark 2.8 (More general unperturbed part). Our model admits an abstract perturbed part 
while the unperturbed operator —A -|- Vq is explicit. The above listed results remain valid 
if we replace —A -|- Vq by a more general operator Cq, as long as it satisfied the following list 
of conditions: 

(i) £0 maps H^fn) to £ 2 ( 1 !) and is self-adjoint on T>(£o) ■= {u & 77^(H) : Bu = 0 on SH} 

(ii) The restriction £o,n of Cq to 77^(0) with boundary condition Bu = 0 on SD n SH and 
Neumann condition on (9n\9n has spectrum cj(£o,n) C {Ao}U [Ai, 00 ), where Aq < Ai, 
Aq is non-degenerate and has a normalized, a.e. positive eigenfunction i/jq satisfying 
7 ^ 0 ,□V’o = AoV'o- 

(hi) Let be the periodic extension of i/’o to H: + k, Xn+i) = , Xn+i) for all 

x' e \Ri,Xn+i S (0, d),/c G T, and xHm on L 2 {JIn)- Let Co^a,N be the 

restriction of Cq with domain (12.51) . Assume that Aq = infcr(£o,n) = inf cr(£o,a,Af) and 
that ~ as well as ■ 

(iv) For any w G H and a n ^ £a,Af(ea;) we have for the spectral infimum the 

bracketing inequality 

^ „min X{C^p k{uj)) 

where 77 ,7 G M, A = 77^ and = 77r n Ta^N- 

(v) With respect to the decomposition © L 2 (n + (A:, 0 )),£o is a block-diagonal operator. 

Remark 2.9 (More general perturbation). Although we have assumed that operators £1 and 
£2 are independent of t, it is possible to treat also the case when these operators depend on 
t, i.e., £1 = £i(t), £ 2 ( 7 )- In this case we should suppose that these operators considered as 
acting from 77^(0) into £ 2 ( 0 ) are bounded uniformly in t. The identity in Assumption lAll 
should hold true uniformly in t, i.e., 

(£i(f)'0o, V’o)L 2 (n) = 0 for each t G [-toTo]- 
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And inequality (I2.7p should be modified as follows: 

(£2(i)i/’o,i/’o)L2(n) “ (t^(0>^iWV’o)L2(n) ^ cq > 0 for each t e [-to,io]- 

where constant cq is independent of t, U{t) is the solution to equation (|2.6p with Ci = Ci{t) 
and U{t) is orthogonal to tpQ in L 2 (n). Then all the above results remain true since their 
proofs remain unchanged. 

The structure of the paper is as follows: The next section presents various specific examples 
which are covered by our general model. They were, in fact, the motivation and origin for the 
choice of the abstract model. Thereafter follows Section 0] with the proof of the variational 
lower bound on the ground state energy on the finite segment and Section [5] with the proof of 
the abstract Combes-Thomas estimate. In Section [6] the proofs of the probabilistic estimates 
are provided. 


3. Examples covered by the general model 

In this section we provide several examples of perturbations covered by our results. Namely, 
we discuss particular cases of operators C{t) satisfying assumptions (lAljl . ()A2p . In what 
follows, we check only this assumptions since they suffice to establish all results presented in 
Section [2j 


3.1. Linear perturbations with positive coupling constants. Condition (Al) imposes 
a quite strict condition on the linear part of the perturbation C. However, if we restrict our 
considerations to non-negative coupling constants then much more general linear perturba¬ 
tion are allowed. To see this, we consider the situation Ci = 0, T 3 = 0. Then (Al) is trivially 
satisfied, (A2) requires cq = {C2'tpo,'ipo) > 0, and C{t) = t^C2- Thus C{£UJk) = £^u:‘f,C2-, hence 
we have non-negative coupling constants as prefactors. Note that any random variable 
Tj: ^ [ 0 , 1 ] can be written as rj = ujq for some random variable oiq: H —)• [— 1 , 1 ], so in 
the case of non-negative random coupling constants the power two is no restriction. In this 
situation Theorem O gives: 

Theorem 3.1. Let 6 G (0, 1), t] € Xfcgr[0, 1], a € T, A^ € M and 

T ^ T 'dkS{k)C 2 S{—k) 

with domain as in (I A 5)) . For sufficiently small 5 this is a selfadjoint operator. 

Then there exist positive constants c\, C2, Ni such that for 

(3.1) > A^i, and 0 < (5 < 

the estimate 

- Ao ^ ^ 

holds true. 

Here A('H^^(r/)) denotes the lowest eigenvalue of The theorem covers the case 

where the random variables are non-negative, the perturbation is linear and in an average 
sense positive. 

In the present situation Corollary 12.71 takes the form of 
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Corollary 3.2. Let a € F, 7 e IN, 7 ^ 17 and ci be as in Ii2. ^) . Choose 5 > 0 and £ M, 
N ^ msix{N^ , KJ , N 2 } , such that N = {5 /. Let /3i,,82, mi,m 2 , Bi, B 2 , C5 he as in 
Theorem 1^.61 Then 


V A ^ Ag + 


1 


7/16 


llXEiiBaNiv) - ^) XB 2 ll ^2 




7/16N 




-n(7-l)/4 


-C4 


Here denotes the distribution measure of the stochastic process rj. 


3.2. Potential. The canonical example is the perturbation by a potential: 

Cit) = tVi+t^V2. 

Here Vi, V 2 are measurable bounded real-valued functions defined on □ and £ 1 , £2 are just 
operators of multiplication by Vi, V 2 - 
Assumption (jAip reads as 

(3.2) J Viipl dx = 0 , 

n 

while Assumption ()A 2 p takes the form 

(3.3) J V2'ipldx> J ViU'ipodx. 

□ □ 

Here Lf solves the boundary value problem described in Assumption (IA2I) with the right hand 
side £i'0o = ViipQ. Since i/’g > 0 is positive everywhere, ()3.2p implies that either Hi = 0 
almost everywhere, or Vi changes sign. It is clear that given Vi obeying (|3.2p . there is a wide 
class of potentials V 2 satisfying (|3.3I) . 


In all the following examples we assume that the boundary condition \2. 31) consists of Dirichlet 
boundary conditions. 

3.3. Magnetic field. The next example is a random magnetic field. The perturbed operator 
reads as 

'H^oj) = {iV + A^f, A^ = e ujkS{k)AS{-k). 

Here A = A{x) = (Ai(a:),... , An+i{x)) is a real-valued magnetic field which is assumed to 
belong to C^PO) and vanishing on the boundary of dDh 

(3.4) A{-,Xn+i) = 0 on dCf for each Xn+i £ (0,(i). 

Since 

(iV + A^f = -A + 2iA^ • V + idiv 
operators £ 1 , £ 2 , £3 are given by the identities 

£1 = 2 iA • V-|-idiv A, £2 = |A|^, £3 = 0 . 
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Let US check Assumption (|Aip . We calculate 

dipo 


(£iV’o,'0o)L2(n) = i j V'o + V'odivA^ dx = i J ^ + ijj^ div dx. 

□ □ 

Now we integrate by parts employing Dirichlet boundary conditions for i/jq and A: 

n 


' I ^"*) ‘‘^='I '>’« 

□ □ 

^ n 

dXn+li’liXn+l) J '^-^{x',Xn+l)dx' = 0 


□' 


since for each Xn+i € (0, d) 

_ r f) A . _ _ /■! /■! 

= 5 ^ / dxi... 

j=i oxj Jo Jo 

which can be checked by integration by parts. 

To check Assumption (|A2h . we first observe that U = iU, where U is orthogonal to ipo in 
^ 2 ( 0 ) and solves the equation 


dAj . 

dXn-;\ - [X ,Xn+l) = 0 

dXj 


Hence, 

(3.5) 

and 


(^n — Ao)U = 2A • Vipo + i/’o div A. 

{U, Cl^po)L2{'0) ={U , 2 A • Vipo + tpo div A)i^^(j2) 

= \\vWL2iu) + {yoU,U)L2in)-^4U\\ Lfn); 

Co = y (|Ap?/)Q — 2C/A • Vi/^o — lAot^div A) dx. 
n 


Let us prove that 
(3.6) 


Co = 


A + V 


U 

TpO 


dx. 


We first observe that functions U and i/jq satisfy the same boundary condition on 50 n SH 
and this is why function ^ is well-defined and belongs at least to Lf^(n). 

To prove (j3.6p . let us calculate the difference 


J := 


A +V 


n 


U 


— |Api/)Q -I- 217A • Vipo + '>PoU div A 1 dx. 


Since 


U U 
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we get: 


J = 


[ ( 2 'il;oA-VU+ \VU\‘^ - 2 ^VU + + '>poUdwA\ dx. 

J \ Wo % ) 


Integrating by parts, we obtain: 


■I 


U 


1 


f2 VV'O 


-2 / —VC/ • VV^o dx = - VU^ ■ —Vipo dx= div dx 


i’o 


V’o 


- AoC/2 - - —ivv^op ] dx. 


□ 


V’O 


Hence, by (13.5p . 

J = j {2^ljoA-dWU+ \VU\'^-VqU^ - AqU^+ 'ilJoUdivA)dx 
□ 

=2 J {i^oA • div C7 + div A + UA-A/i;^)dx = 2 I divipoUAdx = 0, 


□ 


where the latter identity has been obtained by integration by parts. Hence, identity (|3.6p 
holds true and therefore. Assumption (IA2p is satisfied. It means that we can apply the 
results of the present paper to each weak random magnetic held provided (13.41) is satished. 


3.4. Metric perturbation. One more possible example is a random perturbation of metric. 
Here operator reads as 


n+1 


d d 

n%u) = -A + Vb - ^ ^ — {eujkaijix' - k,Xn+i) + e'^ulbij{x' - k,Xn+i)) 


fceP *,i=i 


dxi 


where aij : □ —C, ftjj: □ ^ C, (z, j = 1, ..., n + 1) are complex-valued functions belonging 
to C^pn), vanishing on dlA' x [0, d], and satisfying the symmetry conditions 


aij = aji, bij = bji in □. 


The operators Ci, C 2 , C 3 are given by the identities 

^2 = - ^3 = 0. 


i,j=i “ ‘ hi=i 

Integrating by parts, we rewrite Assumption (lAip as 

n+1 

(3.7) 0 = / 'tpod^i'tpo dx = 


I' 

□ 


n 


d'ipo d'ipo 




/' 


/ dij^o y 
Vdx„+i/ 


dx. 


This identity holds true for a wide class of functions a„+iri+i. The simplest example is 
an+in+i = 0. We stress that ()3.7p makes no restrictions for other coefficients a^j, {i,j) / 
(n + l,n + 1). 
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Assumption 


here looks as 


(3.8) 


^ f dlpo 


bn+l n +1 ^ 


dlpo 

dx. 


n+1 


*J =1 
n+1 

E 

i=l 


dil)o dU 

dxj dxi 


dx 


n+1 


S'i/’o dU 
dxn+i dxi 


dx > 0 , 


where U is orthogonal to + in ^2(0) and solves the equation 


n+1 


{'Ho - Ao)U = - y~] 


dain+1 dipo 


^ dxi dxn+i' 

Inequality (13.8p is satisfied by a wide class of functions +■, Uij. Here the simplest example is 

bn+ln+l ^ 0, CLiji+l — ^n+li — 0 


and other coefficients are arbitrary. (Here we also have to assume that the operator 'Ho does 
not have a constant function as the ground state, as it happens when Vq = 0 and the boundary 
conditions in (12.4p are of Neumann type.) Then the right hand side of the above equation 
for U vanishes and the left hand side in (j3.8l) reduces to f /3n+i n+i ( dx, which is a 

strictly positive integral. 


3.5. Integral operator. The operators £* need not necessarily be differential expressions, 
as above, since we make very weak assumptions in their definition. An example of a non¬ 
differential operator is an integral operator; 

{Ciu){x) = j Ki{x,y)u{y)dy, i = l,2, €3 = 0, 

□ 

where Ki, i = 1,2, are measurable functions dehned on □ x □ and obeying the symmetry 
condition Ki{x,y) = Ki{y,x), (x,y) G □ x □, f = 1,2. 

Assumption (lAll) is equivalent to vanishing of certain mean for Ki: 

/ Ki{x,y)'ipo{x)'tpoiy)dxdy = 0. 

□ xD 


If we suppose Ki{x,y) = K{x)K{y), Assumption (lAlIl becomes equivalent to 


/A-(x+„W<ix = 0 

□ 

and it implies that TiV’o = 0, t/ = 0. Then Assumption (IA2p holds true provided 


/ A' 2 (x, y)'tpo{x)'tpo{y) dx dy > 0 , 

□ xD 

and this inequality is satished by a wide class of kernels K 2 . For instance, the latter inequality 
holds true provided kernel K 2 is non-negative and does not vanish identically. 
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3.6. Boundary deformation. Our next example is devoted to a geometric perturbation. 
Let y = {y',yn+i), y' = (yi,... ,yn) be Cartesian coordinates in IR"’ and iR”"''^, and g = g{y') 
be a non-zero real-valued function defined on IR” and belonging to C'^(R"'). We suppose that 
the support of g is located inside i.e. g vanishes in a vicinity of d\J! and outside □. We 
introduce the function 

9fj{y') = '^s^kg{x' - k,Xn+i). 
fcer 

It is equal to eujkg{y' — k, y-n+i) on k + Employing this function, we define a weak random 
perturbation of the layer 11 : 

:= {y:y'£ R”, g%{y') < yn+i < 9 %{v') + d}. 

The boundary of 11'^ can be regarded as a weak random wiggling of dll. 

In n'" we consider the Dirichlet Laplacian, which we denote by T-L^{ui). The operator TL^{u) 
does not satisfy our assumptions since it is defined on a domain 11 '' depending on a small 
parameter. But it is possible to transform this operator to make it fit our model. Namely, 
one can verify by direct calculation that changing variables x' = y', Xn+i = yn+i — 9 lj{y')i we 
keep the spectrum of 'W{uj) unchanged and we arrive at the operator 

(3.9) =-A-divP^V in R, 

where is (n -t- 1 ) x (n -|- 1 ) matrix defined by 


= 

/.I 


0 


-(V'y5)‘ 


V'-yf 


€ 

|V'yf'2 


^UJ\ 


Vy^ = 


/Ml' 

dx\ 


M 

'dXr, 


The operator (|3.9I) corresponds to (12.2p with 


^. = E 


d dg d d dg d 


+ 


^ dxn+i dxj dxj dxj dxj dxn+i ’ 


€-2 = -\Vgf 


02 


'-'Xn+1 


9 = 9{x'), Ts = 0. 


We proceed to checking Assumptions (lAip . (IA2p . Integrating by parts and taking into 
consideration that g vanishes in the vicinity of OD', we get 


J dx = 


d dg 
dxj dxj 


d'lpodxn+i dx = 0 


n □ 

and Assumption (|A 1 I) is satisfied. 

To check Assumption (IA 2 p . we first observe that 

#0 


d^i'4’0 = 


and, integrating by parts, 
(3.10) 


dXn-\-l 

2 d'^'llJo 


^x'9i 


J '0o^2V’o dx = - y dx = Ao^ 


iV'ypdx'. 


The equation for U reads as 


(3.11) 


{-A-Ao)U = 


d'lpo 

dXn-\-l 


Ax>g 
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and thus, integrating by parts, 


/ 


dXn-\-l 


U A^igdx = 


/ 

□ 

/ 


di>o 

dXfi+i 

dipo 


gA^iU dx 


dx. 


n+l 


- Ao V - -^^x'g 1 dx 


dxl+i 


dx. 


n+l 


f f dipo ^ j f dtpo 

/ \1- gAx'gdx- / - - g 

J VdXn+l/ J dXn+l 


dxl+i 


+ Kn] U dx 


= An / IV 


'g\^dx' - j 


□' n 

Together with (j.S.lOp it implies the formula for cq: 

d^ 


d^ 


dil^o 
dxn+i^ \dxl 


(3.12) 


Co = 


/ 


dipo 

dXn+l V^^n+1 


-9 


‘'n+l 


+ An I 17 dx. 


+ An I 17 dx. 


To check the sign of cq, we solve equation (j3.11l) by separation of variables. Namely, since 


2 . TT 

VO = \/^sm-Xn+i, 

V u d 


we can write the Fourier series 

#0 


dXn+l 


^ ^ lAm(^n+l) •— v d d 


m=l 

Then we represent U as 


(^m. • — 


d'lpo 

dXrL-\-l 


'if^m dx^i-^-X. 


f/(x) — ^ ^ CLfyiUffii^X )'0m(^n+l)? 


m=l 


and obtain that Um should solve the equation 

,.2 

(3.13) 


TT 


-A^> + - l)j Um = A^fg in □' 

subiect to Neumann condition on d\A’. We substitute the above Fourier series for and 

U into (j3.12D to obtain 

2 oo 


(3.14) 


Co = 


vr 


- l)(fi') 17m)L2(n')- 


?n=l 


We represent function 17m as 
(3.15) Um = —I? + Wm, 

and in view of (I3.13F Wm solves the equation 


f-A^-/+ ^(m^ - l)j Wm =-^(m^ - 1)5 in □' 
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subject to Neumann condition on 90'. It yields 


(3.16) 


W/'W iP 




TT 


dP 

TT^ 2 

<^(m 


1)(5'; ^m)l, 2 (n') 

l)ll5'llL2(n')ll^m||L2(n')- 


Hence, 

ll^m||i,2(n') < Ilfi'IUzin')- 

Here we have a strict inequality, since in the case of identity, it follows from (|3.16l) that 
V'Wm = 0 and Wm = const that contradicts equation for Wm- It follows from the obtained 
inequality and (I3.15P that (<7, C4n)i,2(n') < 0 each m ^ 1. Therefore, each term in the 
series in the right hand side of (I3.14p is negative and cq > 0. Thus, our operator satisfies 
Assumptions (HU), (HU and we can apply the results of this paper to a weak random wiggling 
of the boundary. 


3.7. Random operators in mnlti-dimensional spaces. Now we show that our setting 
covers not only operators defined in a hnite-width layer in IR"'^^, but operators defined on 

n 

the whole Euclidean space as well. Recall □' := {x' : x' = ^ Ujej, ai £ (0,1)} C R"'. Let 

i=l 

C'{t) := 

where £( : —)• L 2 ip!) are bounded symmetric linear operators and is bounded 

uniformly in f G [—fo,to]- In n-dimensional Euclidean space we consider the operator 

n'^u:) := -A,, + ^5'(fc)£'(ewfc)5'(-fc), 
fcer 

where A^./ is the Laplacian in R" and S'{k) is a shift operator: {S'{k)u){x') = u{x' + k). 
Then T-L'^ is a random self-adjoint operator in L 2 (R”) with a similar structure as 'H'". The 
only difference is that it acts on functions in R”. 

Based on the £' we define operators £* : R^(n) —)• L 2 (n): 

{Ciu){x' ,Xn+l) = £'u(-,x„+i). 

Thus £( acts on x' i—)• u{x',Xn+i), while x^+i is regarded as a parameter. The result is a 
function depending on x' and x^+i: it is precisely CiU. 

Now that we have Ci, C 2 , and £3 at our disposal, the operator is defined as in (12.21) . 
We choose Neumann condition on SH, Vq = 0, and d = tt. The spectrum of T-L^ can be found 
by separating the variables x' and x„+i. Namely, 

CX) 

(3.17) a{n^)= [J a{n^+ m^) 

m=0 

since we can represent each function u in the domain of T-L^ by its Eourier series: 

00 

u{x) = Um{x') cos mXn+l- 

m=0 

The eigenvalue Aq in our case vanishes: Aq = 0. Assumptions (HU) l!A 2 () take on the form: 
(Al’) The identity f jC'ildx' = 0 holds true, where l(x') = 1 in 
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(A2’) Let U' be the unique solution to the two equations 


-Hu'U' = £il, 



= 0 . 


□' 


Here T-Lu' is the negative Neumann Laplacian on We assume that 


c'o := J C2ldx - > 0. 

□' 

Once these assumptions are satisfied, by (13.1711 and Theorems 12.511^61 we obtain immediately 
the analogues of these theorems for T-L'^. 


Theorem 3.3 (The result described in the introduction). There exist positive constants 
c^, N[ such that for 

(3.18) N > N[, and 0 < e < 

the estimate 

(3-19) - A-o ^ ^ Wfc 


holds true. 


Theorem 3.4. Given 7 G IN, 7 ^ 17, there exist constants C 3 , C 4 , N[ such that for N > N[ 
the interval 


•— 1 ? 4 

_E(|a;fc|)iV2 

is non-empty. For and e ^ In, the estimate 

P (w G H ^ ^ 

holds true, where C 4 depends on pL only. 

In this theorem is the lowest eigenvalue of the operator jy ■ The latter is the 

restriction of to 

n 

Hq,,^ := "[x G : X = o + ^ ^ ciiCi, o-i G (0, iV) ]■ 

i=l 


with Neumann boundary conditions. 


Theorem 3.5. Assume the hypothesis of Theorem \ 3.4\ let e G In and fix / 3 i ,/32 G To^atj 
mi, m 2 > 0 such that Bi := H^^^ C -^2 := m2 ^ noiv T/ien there exists a 

constant C 5 independent of e, a, N, /3i, /32, mi, m 2 such that for N N[ 


UJ G II: \\xBi{H%n - < 2\/iVe' 


Cg dist(S]^ 1-^2) 

Vn 


>l-N 






where || • || stands for the norm of an operator in L2iJl'^N)- 
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4. Deterministic lower bound 

The essential milestone in proving our main result is a deterministic variational estimate 
provided in Theorem 12.II For the reader’s convenience we formulate it here once again. 

Theorem 4.1 (Theorem [2T] above). There exist positive constants ci, c^, such that for 

(4.1) > A^i and 
the estimate 

(4.2) A^^^(a;) — Aq ^ 
holds true. 




C 2 e 




The rest of this section is devoted to the proof of the above theorem. Throughout the proof 
by C we denote various constants independent of e and N. 


4.1. Setup for analytic perturbation theory. We begin with considering operators TL^ jv(^)’ 
i.e., the Schrodinger operator —A + Vq in subject to boundary condition (12.3p on 'ja,N 
and to Neumann condition on 5nQ,^Ar \ 'ya,N- The lowest eigenvalue of operator is jy = Aq 
and the associated eigenfunction normalized in L 2 {Ila.N) is N ~2 2 pQ. Provided N > Ni and 
Ni is large enough, the second eigenvalue is Aq + where x > 0 is the second eigenvalue 
of the negative Neumann Laplacian on Then there exists C > 0 such that the ball 

.B := {A e C : |A-Aol ^ CN-^} 


contains no eigenvalues of TL^ Ar(0) except Aq and the distance from B to all the eigenvalues 
of TL^ j^{0) except Aq is estimated from below by CN~‘^. 

For X £ B \ {Aq} the resolvent {Tl^ — A) ^ is represented as 


(Kn - A)-^ = 


1 (• i'0o)L2(n„,jv) 


Ain 


An — A 


IpO + TZa,N{^), 


where TZa,N is the reduced resolvent. It is an operator from L 2 (nQ,^ 7 v) into Ff^(nQ,,Ar). Its 
range is orthogonal to ^jJQ in L 2 (n.a,N)- Moreover, by analogy with O Lm. 5.2] one can prove 
easily the following lemma. 


Lemma 4.2. For X £ B and f £ L 2 {Ila,N) the estimate 
where C is a constant independent of X, N, f. 


At the next step we describe the minimum of A^ ^ w.r.t. 

Lemma 4.3. (a) For each fixed (sufficiently small) value o/e ^ 0 the minimum of X(^ as 
a function of the variables cok, k £ is achieved for cok = 0, k £ Fq,^ 7 v 
(b) Indeed, there exists p £ (0, oo) independent ofe,a,N, and the configuration D = (a;o)fcgr 
such that 

(4.3) VaGF,AG]N,a;GD,e^0: A^ Ar(<^) ^ min (cqw^ — pjecj^j) 
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(c) Consider the particular configuration w € with = 1 for all k Then 

(4.4) VQ:Sr,A^£]N,e^0: ^a,Af(^) ^ ^0 ~h co£^ — pe^ 

and 

(4.5) Ve ^ 0 : inf a{TC{Cj)) ^ Aq + co£^ — pe^ 
where p is the same constant as in (b). 


Proof. We begin with the case A" = 1. Then ^ is the lowest eigenvalue of operator 
This operator is considered in cell and it given by 


TL%1 — —A + Vq + S{ — 0')£,{£U!a)S{0'). 

By means of regular perturbation theory we can write the first terms of the asymptotics for 




Aa,l =Ao + {C{£UJa)lfo,1po)L2{n) + (^(£Wa)^«,l(Ao)V’O, V'o) 

=Ao + eLOa{C.l'f> 0 ,'lfo)L 2 {U) + {ki-2'f’O “ (Ao)TiV’O; '^o) + C){£^UJ^). 

We apply assumptions (jAlji . (IA2I1 to simplify the above expansion. By (lAljl the next-to- 
leading term vanishes and it is easy to infer form (jA2ji that 

{C2'lp0 — '^l’^o,l(Ao)^lV’ 0 ) = Co > 0 . 

Hence, there exists a constant p independent of e, a.,ijJa such that 

A«,i ^ Ao -h e^ujIcq - p\£^ojI\. 

This identity implies that i achieves its minimum Aq as a function of Ua for iUq = 0. 

We proceed to studying A^ In domain Hq^at we introduce additional Neumann conditions 
on lateral boundaries dOk \ ^H of Dfc for each k £ Ta^N- By the minimax principle it gives 
the lower bound for A^ jy: 

(4.6) A^_jv > min ^ Aq + rnin (cqw^ - p|£w^|) 

feSi a,iV KSia.JV 


At the same time, it is straightforward to check that as Wfc = 0, fe £ Ta,N, the lowest eigenvalue 
of ^ is Aq and the associated eigenfunction is fio. This completes the proof of (a) and (b), 
and (j4.4l) is a special case of (b). For the second bound (j4.5l) in (c) we note that, as above, 
the introduction of additional Neumann boundary conditions yields 

inf a{PP{W))) inf A^ fiu) ^ Aq -|- min (cq - pe). 
fcer ’ fceEa.jv 


□ 


Let us show that A^ belongs to B. In accordance with (14.61) . A^ ^ ^ Aq. By the minimax 


'a,N 

principle we obtain the upper estimate: 
|2 


^a,N ^ 


ll^'^o|li2(n„ jy) + (yofi’o,ipo)L2{D) + {^^fi’o,ipo)L2{n^,N) 


=Ao + 


L2(nQ,jv) 

x; ewfc(TiV'o,V’o)L2(n) + £^((^2 + £Wfc^3(ewfc))V'o, 

k&Ta,N 


L2{^a,N) 
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By assumption (lAip . the sum ^ ea;fe(£i'0o! ^o)L 2 (n) vanishes and we can continue esti- 

k&Ta,N 

mating as follows, 




^a,N ^ ^0 H- 

Thus, as A > A^i and Ni is great enough, jy belongs to B. 


^ An + Ce ^ An + 




4.2. Non-self-adjoint Birman-Schwinger principle. To obtain the desired deterministic 
estimate, we apply the non-self-adjoint modification of Birman-Schwinger principle proposed 
in |10] , in the same way as it was done in m, 13 Sect. 5]. It leads us to the equation for 


(4.7) 


1 


“ Ao — + ^a,NT^a,N{^a,N)) ^^a,Ni’0, i’o) 


where I denotes the identity mapping. As A G B, by Lemma 14.21 and the boundedness of Ci 
we have the estimate: 

(4.8) \\Cl^^na,N{X)\\ ^ CeN\ 

Hereinafter || • || stands for the norm of operators in L 2 {Ila,N)- Thanks to the above estimate 
and (SH), we can conclude that for a properly chosen Ni 

ll^a,iV^a,iv(A)KC<l. 

Hence, operator (X -|- j^TZa,N{X))~^ is well-defined and can be estimated as 

||(X + £^,Aa(A))"'KC 

uniformly in e, N, A. 

4.3. Taylor expansion to third order. Equation (14.71) is the main tool in proving the 
desired estimate for N ~ A-o- We represent (X -|- jv"^oa(A))~^ as 

+ ^%,n'^ol,n{X)) = X — {X) -\- (X^ jv77q,^jv(A)) (X -|- jy77Q,^jv(A)) 

and substitute this representation into (14.7p : 

Xa,N ~ Aq —■^^('^q;,A'V' 0) V’o)L2(nQ,jv) ~ ('^o,IV^a,Af(Act^jv)'^o,Af'0O) V’o)£,2(nQ,jv) 

+ ^ (('^a,7v’^aA(Aa,Ar)) + ^a,N'^a,N{Xa^pq)) ^X^ jyV’O, V’o) 

We rewrite the obtained equation by employing the resolvent identity 

B.a,Af(A“7v) “ '^a,Af(Ao) = {K,N ~ Aq) ^o,Af (Ao)B.a,Af (A"^r) 

as follows: 

Xa^N - Ao =J^{{^a,Ni>0-,'4’o)L2{n^,N) ~ i^a,N'^a,N iX.o)^a,N'^0,'fpo) 12 ( 11 ^,n) 

+ {i^a,N'^a,N{Xa,N)) (^ +(A^^jy)) ^ ^a,N'^0,'^o) 


T 2 (nc,,jv) 


(4.9) 


T (■^a,JV^o,A'(Ao)77Qj7v(A£ jv)'^o,A'V’0) V’o)2^2(11^ jv) 


-1 



22 


BORISOV, GOLOVINA, AND VESELIC 


4.4. Estimates on the individual terms. To estimate the terms in the obtained identity 
we shall make use of the following auxiliary lemma. We recall that Dfc := {x : x — [k, 0) G □}. 

Lemma 4.4. For i G {1,2,3} and each u G the inequalities 


£LOk<S{—k)CiS{k)u, 


^ Ce E lkllLi2(n„ jv). 


kGTn 


Y, eu;kS{-k)CiS{k)i 


fcsr„ 


L 2 (nQ,,jv) 


i' 2 (llc,jv) 

^ Ce ( E Iwfcp ) sup \\u\\H 2 (^Ok)^ 


ker^ 




hold true, where C is a constant independent of e, u, N, u;^. For i = 3, in the above estimate 
we assume C 3 = C 3 {eujk)- 

Proof. Due to the definition of and Cauchy-Schwarz inequality we have 

eujk[s{-k)£iS{k) 


ker. 


Y eojkS{-k)CiS{k)u, 

Y. £^^kiCiS{k)u, 


ker „ 


L2 (IIq,,jv) 


L 2 (n) 


tt. 


fcsr □ 


i 2 (nD 


^ Y ^l‘^fclll^»‘^(^)“lli2(n)ll'^o||i,2(n)- 


Since A: H^iU) L 2 (n) is bounded, we find some constant C such that 
Y ^\^k\\\CiS{k)u\\L^^a)Uo\\L 2 in) ^Ce Y \uJk\\\S{k)u\\jj 2 (□) 


fceEc 


k^Tn 




M E 




kGTn 


Y =<^e( ^kj 

ifesEo, jv / keFc^N 

and we arrive at the first desired estimate. The proof of the other is similar: 

2 


Y^ eoJkS{—k)CiS{k)u 


kev, 


a,N 


i 2 (nc,iv) fcer^, 




Y ^'^^k\\^i^ik)u\\Yn) ^Ce"^ sup \\S{k)u\\jj 2 ^a) Y ^k- 

fcer^.jv fcer„,^ 


□ 


This lemma. Lemma [4.31 and the properties of operators Ci allows us to estimate two terms 
in the right hand side of (|4.9I) . Namely, we have 

1 


m 


(^a,A^a,Af(^o)^Q:,Iv(A"jv)'^a,A'0O, V'o) 


Ce 


€ 


■ keVa^N 

Ce‘^N^ 


Y ^k] '^«,Af(^o)'^a,A(A“7v)T^,wV’0 


H 2 (n„,jv) 


iV" 


fcer^.jv ^ 


(4.10) 
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provided A^i in (14.ip is great enough. In the same way we get 
(4.11) 

\ / L 2 {\-\.a,N) 

= Ce^Ny Y ^iVi E \\m^^^k)Si-k)MUn,) 


kGTot^N 

1 
2 


E ( E £VfcPllV’o|||/ 2 (n) j ^ Y 

' k^T' (y^ N JY JY 


The term iV^ in the third line comes about due to Lemma 14.21 and estimate (14.8D . 
By (|2.1I) we can rewrite two other terms in the right hand side of (14.91) as follows: 


(4.12) V’o)L2(n„,jv) “ ('^a,Af'^a,Ar(Ao)^a,ArV’0, '^o )- Jl + J 2 + h, 


where 


Ji =e E WA:(5(A:)Ti5(-/i:)V'o,'i/’o)i,2(n„.jv) > 

J2=e^ Y ‘^fc(‘5(^)^25(-A:)i/>o,iAo)L2(n„.jv) 

k^T(y^]^i 

-e^ Y Wfca;p(5(/i:)/:l5(-A;)7^„,Ar(Ao)5(p)£l5(-p)'0o,'i/’o)i2(^„^)’ 

p,k&Va,N 

J3=£^ Y ‘^fc(‘5(^)^3(ewfc)5(-A:)i/>o,'iAo)L2(n„.jv) 

k&Ta,N 

E ^fc(‘5(/c)£(ea;fc)5(-/c)7^«,^^(Ao)£^_Jv(^^)^o,V'oj 

A:6r„,jv L2(na,N) 

-e^ E ‘^fc(^a,Ar(^)^a,iv(Ao)5(A:)£(ea;fc)5(-A:)i/>o,i/'oj 

fc6r„,jv 2^i2(na,Ar) 


where £(t) := C 2 + tC 3 {t). 
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Employing Lemmata 14.31 and the properties of /Ij, we estimate J 3 : 


iJsKCe^ uji + Ce^f 






k^Gci 


(4.13) 


+ Ce3f ^ ujU 7^„,A^(Ao) S{k)C{eu:k)S{-h)i^o 

k^GcK^N ksT'(yj\[ 

^Ce^N‘^ ^ ul 

fcSEa JV 




By assumption 


term Ji vanishes: 


(4.14) '^ 1 =^ ^k{CiS{-k)^o,S{-k)^o)L,(u^,^)=^ Y1 = 0. 

fcSEa.JV fcsra,jv 

To estimate J 2 , we shall make use of one more auxiliary lemma. 


Lemma 4.5. The estimate 

(4.15) J 2 ^ e^co ^ ujI 

holds true, where cq is defined in assumption Mgj) . 

Proof. We denote 

/ := ^ ujpS{p)CiS{-p)fio, ux:=TZ{X)f 

P&^a,N 

for A in a small neighborhood of Aq. By ()4.14p . function / is orthogonal to i/’o in A 2 (nQ,^ 7 v)- 
In view of this fact and by the definition of reduced resolvent TZ, it is easy to make sure that 
ux solves the equation 

{Tla - X)ux = /. 

Function ux is orthogonal to fio by the dehnition of TZ. Due to the symmetricity of £1 we 
have 

(4.16) J 2 = e^ ‘^fc(‘5(^)'^25(-/c)V’o,'0o)L2(n,,iv)-e^(«Ao>/)L 2 (^„,Jv)■ 

fcs^c^,iv 

The main idea of this proof is to employ the variational formulation of the boundary value 
problem for ux- Namely, given a domain D and a part 7 of its boundary, by we 

denote the subspace of Lf^(D) formed by functions vanishing on 7 . As A 7 ^ Aq, function ux 
minimizes the functional 

T{u) := f^)L 2 (n„,jv) “ ~ '“)i 2 (n„,A,) 

over and 

-(^^A,/)L2(n,,^) =J^iux)- 

As in the proof of Lemma 14.41 we introduce additional Neumann condition on the lateral 
boundaries dOk \ and it allows to estimate T'iux) from below. More precisely, by W := 

0 iL^(nfc, SDfcnSn) we denote the subspace of L 2 (na,Ar) consisting of the functions such 
k&Ta^N 
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that their restriction on □/£ belongs to hi 911) for each k G It is clear that 

H^{Ila,N,7a,N) C W and thus, 

-iuxJ)L2in^,N) = -^(^a) 


(4.17) 


^STo, w 


The functional in the right hand side of the above inequality is minimized by the solution to 
the equation 

(-A + Vb - Ao) vx = f in Uo,,n \ IJ 90^, 

subject to boundary condition (ESI) on U dDk n 9n and to Neumann condition on 

jv 

U 90^ \ 9n. This solutions reads as vx = t<;fc77fc^i(A)5(A:)£iV'o on Dfc, k G Ta,N- The 

^STq, jv 

restriction of vx on is orthogonal to ipo in L 2 {Ok) for each k G Va^N- Hence, inequality 
(I4.17[) takes the form of 


(4.18) 


- (^A,/)L2(n„.jv) ^ “(^A)/)L2(nc.jv)- 


It is clear that ux and vx are continuous w.r.t. A in L 2 (na,Ar) and = S{k)U on D^, where, 
we remind, function U was introduced in Assumption ()A2p . The latter identity and (I4.18P 
yield 




T 2 (nQ,jv) 


^ col{S{k)U,S{k)CiS{-k)i;o) 


i 2 (nfc) 


= - ^^fc(C'CiV’o)L2(nfe) = “(C'Ci'0o)L2(nfe) ^l- 

fcSPc.JV keVa^N 

We also have 

‘^fc(‘5(^)^25(-A:)V^o,V’o)L2(n„,jv) = Y1 ‘^fc('^25(-A:)V’o,5(-A:)V’o)L2(nfe) 

TV TV 

=(7^2V'o, V’o)L2(n) ^k- 

Now the lemma follows from last two estimates, (I4.16P and Assumption 


□ 

We return back to identity (|4.9I) . First it follows from (|4.10l) and (14.11) that 

1 + ^('^a,Af’^aN(^o)77a,Af(A"Ar)TQ^Ar'0O,V’o)2,2(najv) ^ 2' 

This inequality, (14.lip . (I4.12p . (14.1,ip . and Lemma 14., 5 1 allow us to estimate the right hand side 
of (14.91) and to obtain in this way the estimate for the left hand side: 


^a,N ~ 


y /cgFck^t' 




Ce^N^ E P 

^^Fq. TV 


2 coe 

m 


1 - 


CeN'^ 

2 co 


E 

^^Fq. TV 




By (j4.ip it completes the proof of Theorem 14.II provided ci is small enough. 
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5. COMBES-ThOMAS ESTIMATES 

We establish a Combes-Thomas estimate for the class of operators introduced in Section 
[21 We use that they are block diagonal with respect to the decomposition 0 L2(n + (fc, 0)). 

fcer 

They need not be differential operators. For the reader’s convenience we formulate here 
Theorem 12.61 once again. 

Theorem 5.1 ITheorem 12.61 above.). Let a,/3i,/32 G G, mi,m 2 G IN 6e such that Bi := 
n/3i,mi C B 2 := n^2,m2 C IlQ^Ar. There exists A'"2 G IN such that for N N 2 the bound 

(5.1) IIXB.(WU(w) - A)-AB,llL,(n.,„Hi.(n„„) « 

holds, where Ci, C 2 are positive constants independent of e, a, N, 6, j3i, (32, mi, m 2 , A and 
6 := dist{X, a> 0. 

Proof. We hx w G For arbitrary M we introduce the function J = J{t,M) on [0, M]: 

j{t,M) = t-at)+aM-t), 

where C G C'°°[0,+oo), f{t) = 0 outside [0,1], C^(0) = 1. 

Let ej-, j = 1,... ,n,he the basis in IR” determined by the conditions 

(5.2) {ei,ef)u^=±6ij, 

where 5ij is the Kronecker delta. The sign in the above conditions is chosen by the following 
rule. Given x G Ifa^Ar, we represent x' as 

n 

(5.3) x' = a + 

i=i 

and bj belong to the segments [0, M^W], where Mj > 0 are some constants. The latter 
condition on positivity of Mj determines uniquely the signs in (15.2p and consequently, vectors 
ej-. We also observe that MjN are in fact the distance between the opposite lateral sides of 
the parallelepiped 

By means of expansion ()5.3p and function J we dehne one more function on 11 ^,at: 

n 

Mx') = ^J{hj,MjN). 
i=i 

It is straightforward to check that the gradient of function J* vanishes on the lateral bound¬ 
aries of rtc^Ar and J* G C'°°(nQ^Ar) provided N ^ N 2 and N 2 is large enough. 

Given a > 0, by 7^ we denote the multiplication operator: TaU := e“‘^*u in L2(nQ,^Ar). It 
follows from the aforementioned properties of J* that Ta maps the domain of LL^ ^ onto itself. 
It is also straightforward to check that 

(5.4) r-aAra = A + lpW(a), 

where (a) is a hrst order differential operator whose coefficients are bounded by C(a-l-a^) 
in flc^Ar, where constant C is independent of x G flc^Ar, N, a. 
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The most important ingredient in the proof is obtaining an identity similar to (I5.4p for 
jY- The hrst step follows from the dehnition of 


T-aCl^NTa= 

fcSPct.JV 

Then it is easy to prove the estimates 


(5.6) ||J(- k) J*( ^)||c2(n) ^ T*, 

where constant C is independent of /c € Tq^tv and N. We also mention the inequalities 

e* — 1 ^ fe*, fG[0,+oo), 1 — e* ^ —f, f€(—oo,0 ], 

which can be easily proven by checking the monotonicity of the functions f i-A (f — l)e* + 
1, f € [0,+oo), t ^ t — e“* + 1, f e [0,+oo). These inequalities and (15.61) allow us to 

bound - T(ea;fc)||H2(n„,jv)^i2(n„.jv)- Using the 

expansion (15. 5 p this yield the desired relations: 


||P^^^(a,e,a, A^)||j^2(n^^)^i2(nc,jv) ^ Uae^^, ior V^‘^\a,e,a,N) T-aC%NTa — 


where || • Hx^y indicates the norm of an operator from a Hilbert space X into a Hilbert space 
Y, and C is a constant independent of a, e, a, N. This estimate and ()5.4p imply 


(5.7) T-a'H%]yYa — 'H^j^ + T’{a,s,a,N), 

where P is a bounded operator from H'^{I[a,N) into L 2 (nQ,^Ar) obeying the estimate 

(5-8) ||'P(a,e,a,iV)||/^2(n^ jy)^i2(nc.jv) ^ Uae^^. 

Here C is a constant independent of a, e, a, N. This estimate and the previous identity (|5.7I) 
is the key idea in the proof. It is exactly these two ingredients which allow us to follow now 
the established strategy of the proof of a Combes-Thomas estimate, see e.g. Corollary 3.3 in 

m- 

Our next step is the estimate for the resolvent of We assume that A G [Ao,Ao + 1] 

and provided A belongs to the resolvent set of we have 

ll(7fa,7v - A)"i = ^, 5 := dist(A,a(7f^_^)). 

This identity and the obvious ones 

'^a,N - A = 7Y + i + ^a,N “ i “ A, 

(7fa,Ar — A) ^ = (Tfa^TY + i) ^ + (TIq^at — A — i) ^(7f° 7Y + i) 

yield 

(j 

miN - A)-'llL,(n„,^H772(n„,,) ^ dist(A,u(7f-^^))’ 
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where constant C is independent of e, a, N, and A. This estimate and (15.7p . (15.81) imply the 
inequality 

r'n 

- A)"^|| ^ —-—, 

where constant C is independent of a, e, a, N^ and <5. Hence, for a = C6 with a sufficiently 
small C, 

(5.9) \\['Ha]^ + V{a,s,a,N)—X) || ^ —, 


where constant C is independent of a, e, a, N, and 5. 

Given / 3 i ,/32 G To^at, mi,m 2 > 0 such that C n„,Ar, n^ 2 ,m 2 C na,Ar, by ([531) for 

each normalized vectors G L 2 {Ila n ) we have 

(5.10) 


(lV’l|A:Ri,'71a(^a,iV - ^^A:B2lV'2|)i2(nc,iv) 


^\\T-ainl^ - Xr^Ta 


{T^a,N + ”^(05 ^5 CIj ~ •^) 



where constant C is independent of <5, e, a, N, and A, a is chosen as indicated above, and we 
remind that Bi = B 2 = n^ 2 ,m 2 - Since A is below the spectrum of jy, the integral 

kernel of ("H^ n ~ positive. Without loss of generality we assume that |/? 2 | ^ |/3i|, the 

opposite case is studied in the same way. From now we assume suppi/^j C Bj (j = 1, 2). Then 
it is straightforward to check that 


{\'4’l\XBi,T-ai'Ha,N - ^’^A:R2l^2|)i2(nc,iv) 


^ exp 


1 min J* — max 

V B 2 Si 



• (V’l, ('^a,Af - -^) V2)L2(n„,iv)- 


And since 


minJ* — max J* ^ CdistfHi, H 2 ), 

S 2 Bi 

where G is a positive constant independent of /3i, / 32 , mi, m 2 , a, N, two latter inequalities 
and (I5.10p imply 


I (V^D (^a,7V - V2)L2(n, 




Gi _ 


C2<5dist(Bi,B2) 


where Ci, C 2 are positive constants independent of e, a, N, 6, /3i, / 32 , m-i, m 2 . 


□ 


6. Probabilistic estimates 

In this section we prove our two probabilistic results, Theorem 12.51 and Theorem 12.61 

Proof of Theorem \2.,5[ We follow the main lines of the proof of Theorem 3.1. in [3j. We 
choose A', 7 S M and we let N := K'^. Then up to a set of measure zero we can partition 
Hq^at into smaller pieces 

Ha ,at = 1^ ^/3,K, 

where (J stands for the disjoint union and is the set = KT n Ta,N- We observe 

that the number of elements in the set Ma',^ is equal to {N/K)^ = . On 
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the lateral boundaries of we impose Neumann boundary condition and by the minimax 
principle we obtain 

w ^ min A o . 

l3&MK,-y 

Let us reformulate the above estimate in probabilistic terms. First it implies that 

(6.1) |a; G Q, : “ -^o ^ ^ ^ ^ ^ ~ ^ K 2 1 

/3eMK,-t 


Since random variables cok, k €T are independent and identically distributed, 


(6.2) F(uen: X%K - Ao ^ ^ ')p G L! : A^_^ - Aq ^ iL"2^ . 

l3eMK,^ 

The Cauchy-Schwarz inequality 

fcSPa.K 

and Theorem 14.II yield for K ^ Ni and e ^ ciK~^ 


{w G n : -Ao^iL 2}c^a;Gn: 


(6.3) = <{ w G O : ( Y 

keTa^K 


1 

9\2 ATa 4 
2t2 


y/^e 


— ^ cj G ; 

We choose 6 so that 
(6.4) 


K 


1^,1 K-l 
2^ kfcl ^ 


fcEFo 




. 


C2e 

1^ 


Y 

fcePcK 


1 ^ Art-4 

C <( cu G 0 : —^ > \uJk\^—;=— 

.ir,. 


iF-4 ^ Edcufcl) 
v^e " 2 ^ 


i.e., 

e ^- %■, —^AT-t. 

Y/ctEdwfcl) 

It is clear that this inequality is compatible with (j4.1h provided K ^ ATi, where ATi is large 
enough, depending only on ci,C 2 , 7 ,n, and Edcjol)- We apply the large deviation principle 
analogously as in [U Lm. 4.3]). Hence there exists a constant C 4 > 0 depending on fi only 
such that for each LC G IN 


P 


A k&Fa^K 


E(kol) ) 


< e 


-C4,K^ 


Thus, by (16.3p . (|6.4I) it follows that 


ATT 




keTa^K 


keTa^K 
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as soon as K ^ Ki. Therefore, provided N ^ max{N^, K'^}, where Ni comes from Theo¬ 
rem [4Tl by ()6.2p . (16.Ij) we get 

completing the proof. □ 


Proof of Theorem \2.6i . Here the main ideas are borrowed from the proof of Corollary 3.3 in 
[3]. We first introduce the set 

Htv :=|w e H : - Aq > = jo; G H : dist(Ao, _jv)) > 


= < a; G H : dist(A, a{T-L% p^)) > 


2y/N 


MX G 


Aqj Aq -I- 


2 //vJ j ■ 


For the next step we need a Combes-Thomas estimate as it is given in Section [5l We apply 
it to Hamiltonians j^{uj) with configuration uj in the set Hat. As before we consider 
Bi := c Ila,N, ^2 := n^2,m2 C Hq^at and ^pj G L2(nQ,,Ar) with suppi/'j C Bj {j = 1, 2). 

The Combes-Thomas estimate implies; 

where Ci, C 2 are positive constants independent of e, a, N, 6, /3i, 132, mi, m 2 , and 6 = 
dist{X,a{nij^ioj))),N^N2. 


Now, fix ^ max{Ai^, A^ 2 }) S Hat and A G Aq, Aq -|- 
thus 

KV’I, {T-l-a,N - -^)~V2)L2(n„,jv)| ^ 2Ci'/Ne 2 viv 
By Theorem 12.51 we have bound 


. Then 6 > 


2VN 


and 


nil--] -ciN't 


and therefore, 

PI w G H : VA G 


Aqj Aq -|- 


P(HAr) ^ l-NA^ Pe 


2Vn_ 


|ARi(K,V - A)"^XS2N 2ClViVe -y _ ^n{l l) ^ 


-caNi 


that completes the proof. 


□ 
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